Maximum relative height of elastic interfaces in random media 
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The distribution of the maximal relative height (MRH) of self-afhne one-dimensional elastic in- 
terfaces in a random potential is studied. We analyze the ground state configuration at zero driving 
force, and the critical configuration exactly at the depinning threshold, both for the random-manifold 
and random-periodic universality classes. These configurations are sampled by exact numerical 
methods, and their MRH distributions are compared with those with the same roughness exponent 
and boundary conditions, but produced by independent Fourier modes with normally distributed 
amplitudes. Using Pickands' theorem we derive an exact analytical description for the right tail 
of the latter. After properly rescaling the MRH distributions we find that corrections from the 
Gaussian independent modes approximation are in general small, as previously found for the aver- 
age width distribution of depinning configurations. In the large size limit all corrections are finite 
except for the ground-state in the random-periodic class whose MRH distribution becomes, for pe- 
riodic boundary conditions, indistinguishable from the Airy distribution. We find that the MRH 
distributions are, in general, sensitive to changes of boundary conditions. 
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Recent theoretical progress have highlighted the im- 
portance of extreme value statistics (EVS) in statisti- 
cal physics [TJ [5] . This often yields interesting questions 
of EVS for strongly correlated variables, a field which, 
to a large extent, remains to be explored. For homoge- 
neous systems, it has been recently shown that EVS is 
an interesting tool to characterize the statistics of one- 
dimensional stochastic processes, like Brownian motion 
and its variants, which can be mapped onto models of 
elastic interfaces [3HS]- 

Such elastic interfaces are ubiquitous in nature and 
in the simplest one-dimensional case they can be pa- 
rameterized by a scalar displacement field u(x), for a 
system of length L. Here we consider periodic bound- 
ary conditions (pbc), i.e. u(x = 0) = u(x = L) and 
focus on the displacement, or relative height, h(x) = 

u(x) — (1/L)L da/ u(x'). It can be decomposed into 
Fourier modes 
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where a„'s and 6„'s are random variables, whose statistics 
depend on the precise model under considerations. The 
geometry of such interfaces (IT]) is usually characterized 
by the sample-dependent roughness w 2 and its ensemble 
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where the brackets (...) denotes an average over a„'s and 
& n 's. Of particular interest are self-afnne (or critical) in- 
terfaces, characterized by a roughness exponent £, such 
that Wi ~ I>. Relatively recently, the full distribution 
of w 2 , not only the first moment w 2 , has been consid- 
ered [TT]. It was indeed computed both for Gaussian 
interfaces, i.e. in the case where a„'s and b„'s in Eq. (fTJ) 
are independent Gaussian variables of zero mean and 



variances (a„) = (&„) oc n~ a [as in Eq. (12) below], 
corresponding to a roughness exponent £ = (a — l)/2, 
for a > 1 |llj . as well as for pinned interfaces at the 
depinning threshold [T2J [TJ] , where it was studied both 
numerically and analytically using Functional Renormal- 
ization Group (FRG) [H] . It was subsequently measured 
for a contact line in a wetting experiment, on a disor- 
dered substrate [TS], and also generalized to Gaussian 
"non-Markovian" (a 7^ 2) stochastic processes [TB]. It 
is interesting to note that such Gaussian interfaces can 
be physically regarded as an ensemble of free interfaces 
with non-local harmonic elastic interactions at thermal 
equilibrium, whose Hamiltonian can be elegantly written 
in terms of fractional derivatives of order cv/2 [16 , 



H[{h(x)}] = 
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Here we focus on the maximal relative height (MRH) h n 
of such a fluctuating interface (IT]) defined as 



h m = max h(x) 

0<x<L 



(4) 



This observable was first introduced and studied numer- 
ically for Gaussian interfaces with short range elasticity, 



such that (0%) = (bn) oc n~ 2 [3], i-e. C = 1/2, and it was 
found that (h m ) ~ L^. Then, Majumdar and Comtet 
obtained an exact analytical expression for the full prob- 
ability distribution function (pdf) P(h m \L) of h m and 
showed that it is given by the so-called Airy distribution 

/Airy [HIS], 
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1 



L 1 / 2 



/Airy 



11/2 



(5) 



where /Airy describes the distribution of the area under 
a Brownian excursion on the unit time interval. We re- 
mind that a Brownian excursion is a Brownian motion 
conditioned to start and end in while staying positive 
in the time interval [0, 1]. 

Incidentally, the same Airy distribution appears also 
in various, a priori unrelated problems, in graph theory 
and in computer science 1171 . This result ([5]) was then ex- 
tended to a wide class of one-dimensional solid-on-solid 
models (i , showing the universality of the Airy distri- 
bution ([5]) for interface models with short range elastic- 
ity (and without disorder). The distribution of h m was 
then investigated for Gaussian interfaces with a > 1 in 
Ref. [7j. There it was shown that the distribution has a 
scaling form similar to Eq. ([5]), the scaling variable be- 
ing h m /L^, with ( = (a — lj/2, and where the scaling 
function depends continuously on the parameter a. 

A natural, and very physical, extension of these works 
concerns the maximal height of elastic interfaces in ran- 
dom media: this is the aim of the present paper. In the 
continuum limit the simplest elastic interface model is 
described by the Hamiltonian 



H[{u(x)}] = dx{ C -(Vu(x)f 



V{x,u{x)) - fu(x)j , 



(6) 
equation of motion 
The first term in ([6]) 



together with an overdamped 
dtu(x,t) oc —6H[{u(y)}]/Su(x). 

is the usual harmonic elastic energy, where c is the elas- 
tic constant, and tends to straighten the displacement 
field. The second term is a random potential V(x, u) 
that models randomly distributed defects both in the di- 
rection x along the line and in the direction u along the 
displacements (see Fig. fl for a discretized version of this 
interface to be considered below). The third term repre- 
sents the action of an external uniform field with strength 
/ that tends to drive the interface in direction u. Here we 
consider a random potential which is a Gaussian random 
variable with zero mean V(x, u) = and correlations 



V(x, u)V(x', u') = S(x - x')Rq(u - v!) , 



(7) 



where t-^ means an average over realizations of disor- 
der with correlator Rq(u). We will consider two types 
of disorder which correspond to two distinct universality 
classes: the so-called random-manifold (RM) class, and 
the random-periodic (RP) class. In the RM class Rq(u), 
the bare disorder correlator, is a rapidly decaying func- 
tion over distances larger than a short-distance cutoff a. 



In the RP class, the disorder has periodic correlations in 
the u direction, with a given periodicity ao, such that 
Rq{u + na ) = Rq(u) with n an integer. For a given 
realization of the disorder the steady-state statistics of 
such interfaces is for / = described by a Boltzmann 
weight oc exp [—f3H], with /3 = T^ 1 the inverse temper- 
ature. For / ^ and finite T the system reaches a non- 
equilibrium moving steady-state at large enough times, 
while at T = it can reach a moving steady-state only if 
the driving force is large enough to overcome the barriers 
between the disorder-induced metastable states. 

Such disordered elastic systems ((6| have been widely 
studied during the last twenty years, in particular be- 
cause they have found many experimental realizations, 
ranging from domain walls in ferromagnets [18] . contact 
lines in wetting [151 I19j and fracture experiments [20] . 
They are also relevant to describe periodic structures like 
charge density waves |21j or vortex lattices in type II su- 
perconductors [22J [25]. In these systems, the competition 
between elasticity and disorder leads not only to non- 
trivial ground state configurations but also affects in a 
dramatic way their dynamical properties. In particular, 
when driven by an external force / at zero temperature, 
disorder leads to a depinning transition at a threshold 
value / = / c , below which the interface is immobile, and 
above which steady-state motion sets in. At finite but 
small T, an ultra-slow steady-state creep motion sets in 
below f c , and in particular, at / — 0, the zero velocity 
steady-state coincides with the equilibrium ground-state 
configuration when T — > 0. Interestingly this model m 
has a universal roughness diagram [23] described by the 
crossovers between three types of self-affine interfaces, or 
"reference steady-states" , depending on the amplitude of 
the driving force /: 

(i) at equilibrium (/ = 0) and T = 0, which cor- 
responds to the directed polymer in a random 
medium, where the exponent £ = 2/3 for the RM 
class [25], and C = 1/2 in the RP class [251127] . 

(ii) exactly at the depinning threshold / = f c and T = 
where C ^ 1-25 in the RM class [25] and C = 3/2 
in the RP class [14j EH [29] , 

(iii) in the limit / 3> f c where, in the moving frame, dis- 
order induced fluctuations of the interface become 
an effective annealed noise: for the RM class this 
noise is thermal-like [2U [30] , yielding the Edwards- 
Wilkinson roughness £ = 1/2 corresponding to the 
thermal equilibrium of the system described by 
Eq. (|6| with V(x, u) = 0, while in the RP class 
it is a "washboard" colored noise yielding a larger 
exponent £ = 1.5, identical to the one of the RP 
depinning [26] . l29 ] . 



The third case (iii) effectively describes situations with 
annealed noise rather than quenched noise and in partic- 
ular, for the RM class, the MRH distribution is simply 
given by the Airy distribution (pjl). In the two first cases 
(i) and (ii) sample to sample fluctuations are important 



and one expects different fluctuations of the MRH than 
those produced by annealed noise. It is the purpose of 
this paper to study the MRH distributions in the two 
first cases (i) and (ii) above, that we will denote by RMG, 
RPG for the random-bond (which we also call random- 
manifold) and random-periodic ground states, and RMD, 
RPD for the random-bond and random-periodic depinning 
configurations, respectively. These distributions are thus 
generated by sample to sample fluctuations of the MRH, 
while the ensemble of Gaussian signals (fTl) can be phys- 
ically thought as determined by the Boltzmann weight 
associated to pi with T = 1. 

We show that, in all the situations considered here, 
the distribution of the MRH is very well described by 
the one corresponding to a Gaussian interface, where 
a n 's and 6„'s in Eq. (II) are independent Gaussian vari- 
ables [as in Eq. ( fl2J ) below] of zero mean and variances 
( a n) = (O c* n^^ 2 ^. Despite of this, our numerical 
data show some numerical evidence that these distribu- 
tions are different for the depinning configurations, and 
also for the ground-states in the random-manifold class. 
Note that similar results were obtained for the width dis- 
tribution at the depinning threshold |12j , which were fur- 
ther justified using FRG calculations [13 ]. In that case, 
using FRG close to the upper critical dimension, in di- 
mension d = 4 — e, it was shown that the displacement 
field can be written as u — ^/suq + tu\ where uq is a 
Gaussian random variable of order 0(1) while i*i is a 
random variable with non-Gaussian fluctuations, also of 
order 0(1) [13] . Motivated by these similarities with the 
Gaussian interface, we revisit the analysis of the right tail 
of the pdf of the MRH for Gaussian interfaces and generic 
C where we obtain exact results, using Pickands' theo- 
rem |31j . For the ground-states in the random-periodic 
class with periodic boundary conditions we find that cor- 
rections to the Gaussian independent modes approxima- 
tion vanish at large system sizes and thus the MRH distri- 
bution becomes indistinguishible from the universal Airy 
distribution. Finally, for the equilibrium case we show 
that the MRH distribution is sensitive to the boundary 
conditions, highlighting their importance in the study of 
anomalous, self-affine interfaces. 

The paper is organized as follows: in Section II, we 
remind the results that were obtained for Gaussian in- 
terfaces, and derive the precise asymptotic behavior of 
the right tail of the distribution of the MRH. In Section 
III, we provide the details of our numerical simulations, 
the results of which are then presented in Section IV. In 
Section V, we discuss these results before we conclude in 
Section VI. Some details relative to the use of Pickand's 
theorem have been left in appendix A. 



II. RESULTS FOR GAUSSIAN INTERFACES 

We first start to remind the known results for Gaussian 
interfaces, without disorder. The distribution of h m was 
first studied for the simplest model of an elastic interface 



described by Eq. ((6| with V = 0, / = and periodic 
boundary conditions, i.e. ul — uq. In that case one has 
w| = TL/12. While a first numerical study indicated 
that h m ~ L^, with £ — 1/2, Majumdar and Comtet 
were then able to compute exactly the full distribution 
of h m [HE], P(h m \L) = jjjz / A iry (pt?) where / A i ry is 
the so-called Airy distribution, whose Laplace transform 
is given by 



/Airy(a;)e sx dx 
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fe=i 



2/3o-l/3 
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where a^'s are the amplitudes of the zeros of the Airy 
function on the negative real axis. For instance, a,\ — 
2.3381 . . . , a 2 = 4.0879 . . . , a 3 = 5.5205 ... etc [32]. It is 
also possible to invert this Laplace transform to obtain [U 

w*) = ^fy^^-s/MA&n/* 2 ) , 

X 

n=l 

(9) 
where b n = 2o? r jTl and U(a, b, z) is the confluent hy- 
pergeometric (Tricomi's) function [32] . Its asymptotic 
behavior is given, for small argument by 

/AiryW-Jra^V^- 2 ^ 27 * 2 for x -► . (10) 
The large argument behavior was obtained in Refs. |5"ll3"4"] 



/Airy 0*0 ~72\I~X e 
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In the above expression, the leading Gaussian behavior, 

R 2 , , 

e~ bx , was obtained in Ref. [5] while the sublcading pref- 
actor, 72-t /f x 2 was obtained in Ref. [31] using rather in- 
volved combinatorial techniques [which also allow to ob- 
tain the subleading corrections to ([TT])] . We show below 
that this result (111 can be straightforwardly obtained 



using Pickands' theorem [31j concerning the asymptotic 
properties of the maximum of a stationary Gaussian (see 
also appendix [A] for more details). 

The MRH distribution was then investigated for Gaus- 
sian interfaces displaying a 1// Q power spectrum in 
Ref. [7]. Considering again periodic boundary conditions, 
h(0) — h(L), the relative height field can be written in a 
Fourier expansion as in Eq. (fl]) with a Gaussian proba- 
bility measure: 



V[h] ex exp 



Y^ l 27m ) , 2 , ,2x 



(12) 



corresponding to independent Gaussian random variables 
a„ 's and &„ 's of zero mean and variance, for n > 1 



{a n a n ') = (b n b n >) = 5„ 



21 
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(2irn) a ' 



(13) 



as the one chosen in Ref. [7], which yields in particular 
2L a - 1 



(2tt) q 
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E n ~ a 



(14) 



where C( a ) is the Riemann zeta function. The roughness 
exponent is thus ( = (a — l)/2. In Ref. [7] it was shown 
that the distribution of the MRH takes the scaling form 



P(h m \L) = j-J a . ^ 



(15) 



where from now on, the tilde over any distribution func- 
tion indicates that it corresponds to Gaussian interfaces 
described by Eq. ( 12 ). While the exact expression of the 
scaling function f a (x) is known only for a = 2, see Eqs. 
(p>l [9]). Nevertheless, for a = 2n where n € N, the com- 



putation of the P(h m \L) can be formulated in terms of 
a path integral [7] from which several properties can be 
deduced. In particular, the small argument behavior of 
the scaling function can be studied in detail, yielding [7 
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while the amplitudes B and C can be expressed in terms 
of the smallest eigenvalue eo = £o(a) > of a linear 
eigenvalue problem [7]: 



B = 



1 



2e 



a — 1 



2e 



(18) 



Note that eo(2) = ai/2 1 / 3 so that these formulas (16 
17 18 1 coincide with the small argument behavior of the 
Airy function for a = 2 (10). While this behavior (16) 



where the exponents are given by Eq. (17), was shown 
to hold for a = 2n, numerical simulations demonstrated 
its validity for any a > 1 [7]. We also notice that the 
value of this eigenvalue e (a — 4), corresponding to the 
Random Acceleration Process, was shown to appear in a 
completely different problem of random convex geometry, 
called the Sylvester's question [35] . 

The lar ge argument behavior of the scaling function 
f a (x) ( 15 ) is harder to treat analytically. Numerical sim- 
ulations performed in Ref. [7J showed evidence for the 
asymptotic behavior: 
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As explained in the Appendix [51 we are able to compute 
this exponent 8 as well as the amplitudes D and E us- 
ing Pickands' theorem. One obtains that the exponent 5 
depends on a as follows 



5 = 



-rr ,1 < a < 3 

a— 1 ' 

1 , a > 3 , 



(20) 



while one expects logarithmic corrections for a = 3, which 
thus appears as a threshold value. These values are con- 
sistent with the numerical estimations reported in Ref. 
[7J. They are also consistent with the exact result ob- 
tained in the limit a — > 00, for which 5 = 1 \Q. The 
amplitude E is given by 



E 



(2nT 
4C(a) ' 



(21) 



while the amplitude D can be expressed in terms of a 
Pickands' constant (A3), and has in general a compli- 
cated expression. However for a = 2 and a > 3 it can be 
computed explicitly as 



D = 



72 



C(a-2) (2tt) 



if a = 2 , 

a > 3 . 



(22) 



C(a) 2C(a) ' 



The e xpression of D for generic 1 < a < 3 is left in Eq. 
(A13) in Appendix A. It is straightforward to check that 
the results in Eqs. ( |l9"l [20[ [2l] [22| yield back, for a = 2, 
the result for the Airy distribution ([TT]). 



III. VARIOUS SITUATIONS FOR ELASTIC 

INTERFACES IN DISORDERED MEDIA : 

DETAILS OF SIMULATIONS 

We now describe the various models of elastic inter- 
faces that we study using numerical simulations. Numer- 
ically, it is very convenient to discretize the internal x di- 
rection of the elastic strings described by Eq. Q. Indeed, 
once conveniently discretized, the ground-states and crit- 
ical states at depinning can be obtained by exact and 
fast numerical algorithms. As discussed below, for depin- 
ning we can keep u(x) as a continuous variable. For the 
ground-state it is convenient to also discretize u{x) and 
work on a lattice. None of these practical choices change 
the universality classes we analyze, and at large enough 
L the numerical results for the lattice model should be 
indistinguishable from those of a continuum model. 



A. The ground state 

When the external force is zero and the finite size 
pinned interface is allowed to equilibrate at zero tem- 
perature it reaches the ground state configuration, found 
by minimizing the energy J6J) for a given sample of dis- 
order. Finding such a state for each disorder realization 
poses a global optimization problem which can be solved 
by an exact transfer-matrix method. This method allows 
to find the ground-state of a line on a lattice where both 
u(x) and x are discrete variables. For simplicity we im- 
pose a hard constraint on the local allowed elongations, 
\u(x+l) — u(x)\ < 1. This simplification does not change 
the equilibrium universality class which is thus only de- 
termined by the disorder correlator. This is at variance 
with depinning configurations, where such constraints, or 




(4,6) 



FIG. 1. The plain blue line represent a path on the lattice 
when starting in jo = 3. The pointed maximum is reached 
at i — 4 and its value is ji = 6. If (ji)o<i<L is the sequence 
of heights visited by the path, one must substract the spatial 
average (u) = ^nJi/L to obtain the MRH of this sample. 



more general non-linear corrections to the standard har- 
monic elasticity, can change the depinning universality 
classes [53J |3S] . We thus simulate the directed polymer 
in a random medium as a simple model for the continuum 
interfaces described by J6J), with the universality classes 
determined by the nature of the discretized disorder po- 
tential correlations. 

We work on a finite discrete space: we use a rectangu- 
lar lattice of size L x £. Points of the lattice are indexed 
by the couple of integers (i,j) in the ranges < i < L 
and < j < t. The disordered potential is now a random 
variable Vij given at each point (i, j) of the described lat- 
tice, with zero mean 



Va = o, 



and correlations 



VijVu = Si k Sji 



(23a) 



(23b) 



In our simulations, we choose Gaussian random variables, 
although the precise distribution should not be relevant, 
provided that it remains narrow. We work with periodic 
boundary conditions both in the u and in the x direction: 
the sites (0, j) are identified with the sites (L,j), for all 
< 3 ' < I, and in the u direction as well: (i, 0) = (i,£), 
for all < i < L. A path C of length L is given by the 
sequence of visited sites (i,j), for < i < L— 1, and due 
to periodic boundary conditions, it can be regarded as a 
close path on a torus. Its energy is the sum of disordered 
energies at each visited site: 



E(C)= £ Va 



(24) 



(i,j)ee 



At each step, the path can go from (i,j) to one of the 
three possibilities (i+l,j — 1), (i + l,j) or (i+l,j+l): we 



then select the optimal path amongst the ensemble Cj a of 
paths starting in (0, jo) and ending in (L,jo) = (0, jo) by 
a recursive method [25]. Calling C* this path, one has 



E opt (j )=E(C* )= mm E(C jo ) 

C 30 



(25) 



the energy for the optimum path, minimizing the total 
energy for a fixed starting point jo- However considering 
only these paths C* o , even averaging over the disorder, 
yields a non-stationary two-point correlation function of 
the displacement field. One obtains a stationary correla- 
tion function by considering the optimal path C* among 
all closed paths on the torus. Then 



E GS = mmE(C) = E{C*). 



(26) 



From our ensemble of paths C* , one has to minimize by 
varying the starting point jo on the u axis: 



E GS = min E opt (Jo) 



(27) 



Recalling that L and £ are the lengths of the disordered 
substrate in the x and u directions respectively, one has 
to choose £ > L^ GS in order to explore a region of corre- 
lated paths. In this case, £gs stands for the ground state 
roughness exponent, which can take the values £rmg = 2/3 
[25] for the random-mainfold ground state and £rpg = 1/2 
for the random-periodic ground state, respectively. In 
our simulations, we have taken £ = L to produce ground 
states in the RM universality class, while we took £ -C L, 
for ground states in the RP universality class. The last 
choice is equivalent to construct a quenched potential 
with periodic correlations in the direction of u with pe- 
riod £. In particular, if we take / C i but still large 
compared to the unity, £ controls the crossover between 
the RM and RP universality classes, from £ RM g to £rpg- 
For this reason, in most of our simulations for the RP 
class, we set I ~ 0(1), thus minimizing the crossover or 
finite-size effect. 

We have produced a number of 10 6 samples of each 
lengths L = 128,256,512,1024,2048,4096 in order to 
construct the histograms of the maximal relative height 
h m . The results we obtain for the MRH and discuss in 
the next section are computed from stationary paths C*, 
and are referred with the index RMG, RPG for the random 
manifold and random periodic universality classes. We 
discuss the results for the ensemble given by C* (with 
the arbitrary choice jo = 0) having a non-stationary cor- 
relation function in the section |IVD| for the random- 
manifold case using the index RMG' to distinguish from 
its stationary counterpart RMG. Both cases, RMG and RMG' , 
have the same roughness exponent £ RMG = Crmg' = 2/3. 



B. The depinning transition 

The critical configuration u c (x) of an elastic interface 
at the depinning transition is an extreme solution of the 



overdamped equation of motion 



du(x,t) 
dt 



SH[Mx)}] 
Su(x) 



+ / = o, 



(28) 



such that u — u c for the largest possible force / = f c 
satisfying the above equation, with f c the so-called crit- 
ical force. Above f c , the sum of all forces in the second 



term of Eq. ( 28 ) is always positive, and thus a zero ve- 
locity steady-state solution can not exist. Middleton's 
theorems 128 assure that such a solution exists and it is 



unique. In order to solve numerically Eq. (281 for a given 



realization of the random potential we can discretize the 
string in the direction x into L pieces, thus converting 
Eq. ( 28 1 into an inversion problem. Solving the resulting 



L-dimensional system of Eq. ( 28 ) for large L by standard 



general methods is however a formidable task, due to the 
non-linearity of the pinning force. On the other hand, 
solving the long-time steady-state dynamics at different 
driving forces / both below and above f c is very inef- 
ficient due to the critical slowing down near f c . Fortu- 
nately, this problem has a particular analytical structure 
that allows to devise a precise and very efficient algo- 
rithm allowing to obtain iteratively the critical force f c 
and the critical configuration u c (x) for each independent 
disorder realization without solving the actual dynamics 
nor directly inverting the discretized version of Eq. ( 28 ) 



[2"5l 156] . In this paper we use such method with periodic 
boundary conditions in all directions. This guarantees 
that the critical configurations have spatially stationary 
correlation functions. 

We implement the algorithm to find the critical con- 
figuration in 1+1 dimensions as in Ref. [28 . We dis- 
cretize the space in the x direction while keeping u as 
a continuous variable. A potential V(x, u) satisfying 
d7|) is modeled with random cubic splines. We con- 
sider periodic boundary conditions in both directions, in 
a system of size L x I. When I is large enough, the 
critical configuration sample averaged width is well de- 
scribed by w L = [u c 2 — u^ 2 ] 1 / 2 ~ £,£ with C = Crmd ~ 
1.25 |28| the random- manifold depinning roughness ex- 
ponent. We will denote this Random-Manifold depin- 
ning case as RMD. For small I the average width is well 
described by the random-periodic depinning roughness 
exponent C = Crpd ~ 1.5 [H [Ml US]- Accordingly, 
we will denote this Random-Manifold depinning case as 
RPD. More precisely, for this system it was found that 
w L ~ G(£/L c ™)L Crhd for all values of t with G(x) a non- 
monotonic function of x, such that G(x) ~ 2;( 1_ Ci™/CRm) 
for i C 1 and a with a very slow, approximately loga- 
rithmic, growth for x 3> 1 |29j . Changing the transverse 
periodicity thus allows to crossover from the random- 
manifold to the random-periodic universality class. By 
using this method we can sample critical configurations 
belonging to these two classes, RMD and RPD, respectively, 
and tune the non-universal prefactors w^/L^™, w^/L^" 
to different values. Once the critical configuration u c (x) 
is obtained for each case, we substract the center of mass 
position h{x) — u c {x) — (1/L) L dx' u c {x') and calculate 



the MRH from Eq. Q . Repeating this procedure for dif- 
ferent disorder realizations gives access to the MRH prob- 
ability distribution. In our simulations we use between 
10 4 and 10 5 critical samples of sizes L — 128, 256, 512. 



IV. NUMERICAL RESULTS 

We present numerical results of the probability dis- 
tribution function Pk(h m \L) of the maximal relative 
height h m where the subindex k = RMG,RPG,RMD, RPD 
denotes the different cases: random-manifold and 
random-periodic ground states, and random-manifold 
and random-periodic critical configurations, respectively. 



A. Scaling of h m with L 




FIG. 2. For the random-manifold and random-periodic 
ground states, k = RHG, RPG, and critical configurations 
k = RMD, RPD, we plot the rescaled distributions fk(z) = 
L < " k Pk{z x L^ k \L), obtained from the histograms Pk(h m \L) 
computed for different system size L as indicated in the key. 



For all cases, our numerical data are compatible with 
the scaling law 



P k (h m \L) = -L f k 



fl 7n 



(29) 



where Cfc is the roughness exponent of the case k, Crmg = 
2/3, Crpg = 1/2, Crmd » 1-25, and Crpd = 3/2. In the 
continuum limit L ^> 1, the rescaled functions fk are 
expected to depend only on the rescaled maximal rel- 
ative height z = h m /L'' k . In Fig. one observes that 
the rescaled distributions for different L of the case k in- 
deed collapse for different sizes of the same class, strongly 
supporting this scaling relation (29 1. The different cases, 
RPG, RMG, RMD and RPD in order of roughness exponent, 



have very different scaling functions fa in this L^ scaling. 
In particular, it is clearly visible that the most probable 
value, the mean and also the standard deviation of each 
fk decreases when the roughness exponent increases. Be- 
sides these facts, the curves are somewhat similar, and 
we are interested in testing other simple scaling forms to 
compare these distributions Pk(h m \L) on the same foot- 
ing. 



B. Average and sigma scaling 
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FIG. 3. For each case k = RPG, RMG, RMD, RPD we show the 
average-rescaled distributions <t>k(x) — h m Pk(xh m \L). Differ- 
ent system lengths are shown, as quoted in the key, and the 
rescaled functions collapse onto the same curves, having unit 
average. 



Two other scaling forms have been used in the litera- 
ture [7J [T2] , namely the average-scaling and the sigma- 
scaling. The average scaling is defined by 



Pk(h m \L) 



<i>k 



where the average is computed as 



lim — 



dh m h m Pk(h m \L). 



(30) 



(31) 



The rescaled variable x = h m /h m has by definition a 
unit mean. It differs from the precedent scaling by the 
non-universal prefactor A^ in h m = Aj c L^ k . This pref- 
actor can be indeed physically relevant if there exists 
a crossover to the asymptotic roughness regime which 
characterizes the universality class. In such cases, if the 
crossover to the asymptotic regime takes place at a char- 
acteristic length L x , we can expect h m « A x (L/L x ) l ' k , 
with A x ~ h m (L x ) a characteristic maximal height for 



10" 



10" 



10' 



10' 



10"' 



10" : 









" 


• 






- 


A W RPG (y) (L=4096) 

" i " VrmgM (L=4096) 

i *BMD W < L =256) 

i* l'RPD(y) (L=256) 

- J ^Airy(y) 


• 

* 
■ 




_ 


j" 









-2 



1 2 

y=( h m - R m) /o 



FIG. 4. For each case k = RPG, RMG, RMD, RPD wc show the 
sigma rescaled distributions ipkiy) = oPk{h m + cry\L), ob- 
tained with L — 256 for the depinning cases and L = 4096 
for the ground state cases. Other curves for different L col- 
lapse on theses ones and are not shown for clarity. The Airy 
distribution is shown also with sigma scaling. 



a system of size L x . We thus obtain A^ w A x /L x . 
Such crossovers usually can depend on microscopic de- 
tails, such as the strength of the disorder and the elastic- 
ity (e.g. the Larkin length), temperature or the spatial 
discretization when the interface is defined on a lattice. 
In the case of moving interfaces such crossovers can in 
addition depend on the velocity [23]. As a peculiar case, 
in very elongated samples the prefactor Ak could also 
depend on the transverse dimension of the system if the 
configuration carries with it an extreme value over the 
sample, such as the depinning threshold for the critical 
configuration [35]. We present this average scaling in 
Fig. [3] We see again a very good collapse for the differ- 
ent system sizes within each case k, and also observe that 
the scaling functions 4>k look still very different because 
the width of the rescaled distributions <f>k appreciably in- 
creases when the roughness exponent increases. 

Yet another way to present the data is the sigma- 
scaling defined by 



Pk(h m \L) = -%l> k 



(32) 



where the standard deviation a is obtained when averag- 
ing over the disorder 



(33) 



The rescaled variable y — (h m — h m )/a has zero average 
and unit standard deviation. This is the most general lin- 
ear transformation that can get rid of model dependent 
amplitudes, such as the aforementioned A^. It allows 
then to compare the shape of the distributions on the 



same footing. Actually, the differences between the dis- 
tributions tpk are minute, and not visible on a plot with 
linear axis. In Fig. [4] we show the sigma-rescaled dis- 
tributions ipk obtained from samples with L — 4096 and 
L = 256 for ground states and depinning respectively, us- 
ing log-linear axis (distributions obtained from other sys- 
tem sizes collapse on the same curves and are not shown 
for clarity). One can observe that the ipkS are eventu- 
ally different: this is especially visible in the left tail. 
It indicates that the rescaled distributions remain sensi- 
tive to the roughness exponent £. The Airy distribution 
function (also plotted in sigma scaling) seems to collapse 
very well with the RPG case (see the discussion in para- 
graph IV C below) . Moreover it gives an approximate 
description of the remaining cases RMG, RMD, RPD, furnish- 
ing an idea of the small and large argument behaviors. 
In particular we have checked that the large argument 
behavior of ipk{y), for each case k = RPG, RMG, RMD, RPD, 
is well fitted by a Gaussian tail, ipk(y) oc exp(— d^y 2 ), 
although a precise determination of the coefficients dk 
remains a hard task numerically. Considering this ob- 
servation together with our exact analytical results for 
Gaussian interfaces obtained above ( fl9| ), we conjecture 
that ipk(y) °c ex p( — dky 2 ) is the exact leading behavior 
of ipk(y) for large y. 
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FIG. 6. Right tails of the MRH distributions in average scal- 
ing. We show the data for Gaussian interfaces (open signs) 
with corresponding values a = 2,7/3,7/2,4, and their exact 
asymptotics (continuous lines) computed from Eq. (19 1. We 
also show the disordered interfaces (filled signs) with (t for 
the four studied cases with k = RPG, RMG, RMD, RPD. For each 
case the fitted asymptotic tails as exp(— e^x 2 + c^) are also 
presented (discontinuous lines). 



C. Comparison to Gaussian Independent Modes 
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FIG. 5. Comparison between the distributions <j>k{x) (sym- 
bols) of the MRH of the elastic interfaces in random media, 
for each case C, k - 1/2, 2/3, 1.25, 1.5 for k = RPG, RMG, RMD, RPD, 
obtained from the length L = 4096, 4096, 256, 256, and the 
distributions 4>a {x) (lines) of the MRH of Gaussian interfaces 
with a = 2, 7/3, 7/2, 4 and corresponding roughness expo- 
nents C,k- All distributions are shown with average scaling. 
The Airy distribution, corresponding to C — 1/2, is plotted 
as a reference. 

A natural question arises when studying disordered 



systems: is there a simpler model that can absorb the 
complexity induced by the spatial randomness? If a sat- 
isfying answer comes out, it will permit to obtain more 
information on the disordered system. Indeed, models 
without disorder are often easier to work on numerically, 
and are also more tractable analytically. For our prob- 
lem a good candidate is the Gaussian independent mode 
interface as it can easily describe the self-affine geom- 
etry of the pinned interface. Physically, the sample to 
sample MRH fluctuations due to different disorder real- 
izations are thus approximated by the "thermal" MRH 
fluctuations of a free elastic interface with a non-local 
elasticity as described by Eq. ^. Comparing with Gaus- 
sian signals hence allows to separate the purely geometri- 
cal features of self-affine disordered interfaces, which can 
be successfully described by independent Fourier modes, 
from the specific non-Gaussian corrections generally ex- 
pected from the complex interplay between disorder and 
elasticity. This analysis is thus experimentally relevant 
as it can provide very specific information from a purely 
geometrical analysis of interfaces. 

Hence we compare our results to the distribution of 
MRH obtained with Gaussian interfaces generated via 
the Fourier expansion (nj) with the probability measure 
given in ( 12 ). We therefore adjust the parameter a which 



parameterizes the Gaussian signal to get the correspond- 
ing roughness exponent. Since the roughness exponent of 
the Gaussian interface is £ = (a — 1)/2 we will use the no- 
tation a RPG = 2, a mG = 7/3, a RM D = 7/2, and a RPD = 4 to 
characterize the Gaussian cases used to compare with the 
corresponding disordered interfaces. Average and sigma 



scaling of the MRH distribution is used to adjust the am- 
plitude of the Gaussian interface (or "temperature" ) and 
get the best Gaussian approximation for the ensemble of 
pinned interfaces. 

To start the comparison, we have computed the MRH 
distribution (j)^ in average scaling, for the Gaussian sig- 
nals corresponding to the four cases. In Fig. [5] one can 
see that the distributions for the elastic interfaces in dis- 
ordered media are well approximated by their pure Gaus- 
sian counterpart, especially for k = RPG,RMG. In Fig. |6j 
we zoom on the large x behavior in a log-linear axis, and 
plot the exact asymptotics for Gaussian interfaces which 
we computed above in Eqs. (19 20 21 22 1 - up to a 
rescaling of h m - and the fitted asymptotics for the dis- 
ordered interfaces. As mentioned in the analysis in sigma 
scaling of Fig. [4j the data shows a good agreement with a 
Gaussian tail exp(— e^x 2 ), at least at leading order. One 
sees that the coefficients e& are slightly different from 
those entering the exact asymptotics of Gaussian inter- 
faces. We also observe that the value of e& is closer to its 
Gaussian counterpart for the case of ground-states, RPG, 
RMG. 

To go further and characterize better the possible 
differences with Gaussian signals we proceed as in 
Ref. [12] for the width distribution and compute (numer- 
ically) the cumulative MRH distributions for each k — 
RPG, RMG, RMD,RPD for different sizes L and their corre- 
sponding Gaussian counterparts (denoted with the tilde) 
as references, computed with the largest size L = 4096 
to minimize finite-size effects: 



F k (x) = / dx' M*'), 



(34) 



F k (x) = / dx' cf> ak (x'), L = 4096, (35) 

Jo 

and analyze the difference AFk(x) = Fk(x) — F k (x) for 
different values of the length L. In Fig. [7] we can see 
that for the depinning cases, RMD and RPD, the differences 
AFk(x) seem to saturate for all values of the rescaled 
variable x = h m /h m . This is consistent with what was 
found for the same difference regarding the width distri- 
bution for critical configurations at depinning |12j . In 
Fig. [8] we show the ground-states cases. While the RMG 
case maintains a finite difference as L increases for the 
RPG the differences decrease towards zero for all x with 
increasing size. 

To quantify the difference between the MRH distri- 
butions of elastic interfaces in random media and their 
corresponding Gaussian signal we start with a very sim- 
ple statistical test. We assume that disordered elastic 
interfaces can be mapped onto a Gaussian signal by a 
simple rescaling of its amplitude by a factor \i 



h{x) = x h ( x )y 



(36) 



and thus the two ensembles would be equivalent. This is 
indeed motivated by the analytical predictions for critical 
interfaces at depinning showing that the displacement 



0.04 
0.03 
0.02 
0.01 



-0.02 







AF RMD (x)(L=128) 








A f rmdM (L=256) 








AF RMD (x)(L=512) 


" 


* - 




AF RPD (x){L=128) 




: 




A F RPD (x) (L=256) 




•a °rf. 

fa "■• 




AF RPD (x)(L=512) 




•a °* 
-a a** *Jl 


















- 

















FIG. 7. Differences Aft(i) of cumulative distributions be- 
tween disordered interfaces and Gaussian interfaces, for the 
depinning cases k = RMD, RPD, for system lengths L = 
128,256,512 for the disordered interfaces data, in average 
scaling. 
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FIG. 8. Differences AFk(x) of cumulative distributions be- 
tween disordered interfaces and Gaussian interfaces, for the 
ground-state cases k — RPG, RMG, for system lengths L = 
1024, 2048, 4096 for the disordered interface data, in average 
scaling. 



field can be written as u = \/euo + eui where uq is a 
Gaussian random variable of order 0(1) while ui is a 
random variable with non-Gaussian fluctuations, also of 
order 0(1) [12]. Since h(x) and h(x) describe a self-affine 
interface with the same roughness exponent £ we must 
have, in particular, 



hm = aL c , (h m ) = 5i c 



(37) 



\J(h m - h m y = bL<, ^{{h m -{h m )Y) = bL<. (38) 
For the particular case RPG the parameters of the Gaus- 
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sian signal a and b corresponding to M can be obtained 
analytically, since in this case the MRH distributi on is 
Airy distributed, since £rpg = 1/2. We get a RPG = y/n/8, 
6rpg = ^/5/12 — tt/8. For the other cases, the Gaussian 
parameters can be obtained numerically for a finite num- 
ber of modes. In order to avoid mixing size effects present 
in both, the disordered interfaces and Gaussian signals 
through the number of modes, we have evaluated a and 
b for a very large number of modes. We find indeed that 
the value a/b converges faster to its assymptotic value for 
larger roughness exponents. We have thus fixed L = 4096 
modes for the Gaussian signals, assuring that the values 
a and b are almost numerically converged for the C, = 2/3 
case (and thus for all higher exponents). By making the 
hypothesis (|36|) we get 



h m = x(hm) =>■ X = ~ 
a 



(39) 
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FIG. 9. Values of ratios A^ for each case k = 
RMG, RMD, RPG, RPD plotted with increasing size L for the dis- 
ordered interface data. The inset presents the power law be- 
havior of |Ar PG | towards zero, in log-log scale. 



implying that b/b = a/a if the hypothesis is true. To 
quantify the possible differences we can thus define the 
ratio 



(a/a) ~ (b/b) 
min I a J a, b/b ) 



(41) 



for each case k. In Fig. [9] we show the evolution of 
Afe with the system size L for all the cases, RMG, RPG, 
RMD, and RPD. As discussed above finite size effects come 
from a and b and not from a and b. We can see that 
Afc < for the random-periodic cases (k = RPG, RPD) 
and the ground-state of the random-manifold (k = RMG), 
while Afc > for the random-manifold at depinning 
(k = RMD). It indicates that the variance of the average- 
rescaled MRH distribution for the disordered interface 
data k = RPG, RMG, RPD is larger than the variance of 
the corresponding average-rescaled MRH distribution for 
the Gaussian data, and the opposite is true for the 
RMD case. Focusing on the size-dependence, we see that 
Armd quickly saturate to values of order of 0.03 — 0.04, 
while both |A RMG | and |A RPD | slightly increase and should 
converge as well to a non-zero value for large systems 
L — > oo. On the other hand |A RPG | slowly decreases 
with L and show no tendency to saturation towards a 
finite value. Interestingly, as we show in the inset of 
Fig. [9] we have that for the RPG case |A RPG | w L^ 0Rpa , with 
8 RP q sa 0.5. We find that this power law behaviour is not 
proper to the disordered interface however: the corre- 
sponding Gaussian signal also follows a very close power 
law as a function of the number of modes, as its MRH 
distribution approaches the assymptotic Airy distribu- 
tion. To better quantify the non-Gaussian corrections we 
have performed a two sample Kolmogorov-Smirnov sta- 
tistical test of equivalence between the MRH of sampled 
ground-states and their Gaussian approximations for the 
RMG and RPG cases. We have used the data sampled at 



the largest system sizes L = 2048, 4096 for both the dis- 
ordered interfaces and Gaussian signals in order to re- 
duce at maximum the possible finite-size effects in both 
processes. For the RMG we find a statistically significant 
difference between the distributions, with a probability 
less than p < 10 -8 that the distributions are the same 
(null hypothesis). For RPG however, we find a probability 
or significance of order p = 0.5 for L = 2048, and p = 0.8 
for L = 4096. It is thus very probable that the two distri- 
butions are the same one. Since the Gaussian case tends 
to the Airy distribution as L — > 00 we conclude that the 
MRH statistics in the Random-Periodic class is in the 
same class as the one of periodic normal random-walks 
with C = 1/2. In this respect we note that this extends 
the study made in [B] made for a family of thermally 
equilibrated Solid-on-Solid models without disorder, to a 
disordered system at T = 0. 



D. Boundary conditions 

The dependence on boundary conditions is a relevant 
issue. On the one hand it is experimentally usually dif- 
ficult to realize prescribed boundary conditions for self- 
affine interfaces. In this respect a theory describing "win- 
dow boundary conditions" have been recently proposed 
in which the width statistics of segments become univer- 
sal and independent of boundary conditions |16) . On the 
other hand, boundary conditions usually add technical 
difficulties to the analytical approaches. Numerical stud- 
ies with different boundary conditions are thus interest- 
ing, as they allow us to know the sensitivity to boundary 
conditions of different observables. Here we will address 
their effects on the MRH statistics for particular cases. 

So far we have only analyzed self-affine periodic sig- 
nals h(0) = h(L) with a given period L, either origi- 
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nated from independent Fourier modes or from the inter- 
play between local elasticity, disorder and a driving force, 
which all have stationary spatial correlation functions. 
This means for instance that for any two observables 
0i(x) and 2 (x) the Gaussian signals considered satisfy 
(Oi{x)0 2 (x')) = (Oi(0)O 2 (x'-x)), while the considered 
critical configurations at depinning and ground-state con- 
figurations satisfy 1 {x)0 2 {x') = O 1 (0)O 2 {x' - x). In 
particular, this implies that the MRH can occur, with 
uniform probability, at any position x. 

Non-stationary spatial correlation functions can be 
easily generated using independent Fourier modes, while 
keeping the self-affine geometry, by applying different 
boundary conditions. For a self-affine signal h(x) with 
roughness exponent £ we can for instance impose fixed- 
ends h(0) — h(L) = by constructing a pure sine series 
h(x) = J2 n =i c n sm(nirx/L) with normally distributed 
uncorrelated amplitudes c n — N[0,(2/irn)(L/irn)''] [37]. 
It is easy to see that these signals have in general non- 
stationary spatial correlation functions. Consider, for in- 
stance, the observables Oi(x) — 2 (x) — h(x) — h(0): the 
correlation function (Oi(x)0 2 (x')} vanishes for x = (or 
x' = 0) and x — L (or x' = L). In particular, it is clear 
that for fixed-ends the location of the MRH of each sig- 
nal is not uniformly distributed along L: instead its pdf 
is peaked around x = L/2. The MRH distribution of 
these kind of non-stationary signals is not known exactly 
in general. 

Ground-state configurations of elastic interfaces in ran- 
dom media can also have non-stationary spatial correla- 
tions functions in presence of certain constraints. Indeed 
we have seen that an optimal ground-state path with 
h(0) = h(L) [note that we do not consider in this paper 
"tilted" signals with /i(0) ^ h(L)] can be generated in 
two ways. The one we have analyzed so far comes from 
a minimization over all possible paths in the disordered 
substrate with h(0) = h(L), and periodic boundary con- 
ditions in the displacements direction. This yields sta- 
tionary ground-state configurations. On the other hand, 
if we do not seek the minimum among all possible paths 
such that h(0) — h(L), but fix the extremity h(0) — h(L) 
to a particular value we obtain ground-state configura- 
tions with non-stationary spatial correlation functions. 
In the two cases however, the self-affine structure is pre- 
served, as it only depends on the variance of Fourier 
modes amplitudes rather than in its phases. Finally, let 
us note that the critical configurations at depinning we 
have analyzed have, by construction, stationary spatial 
correlation functions. 

We will address here the effects of boundary conditions 
in the MRH distribution. To this purpose we focus on 
the well known directed polymer in a random medium, 
equivalent to the RMG case for our disordered elastic man- 
ifold. To illustrate the effects in this case we generate two 
particular types of boundary conditions by using the two 
methods outlined above to obtain both stationary and 
non-stationary ground-states with Q = 2/3. One may 
write the height field as a fully periodic series, contain- 



ing both cosines and sines as in Eq. (TTT) in the stationary 
case and as a sine series in the non-stationary case. We 
note that these boundaries conditions would give equiva- 
lent results for £ = 1/2 37 , which is the case of the RPG 
disordered interface. 

In Fig. [lO] we plot the difference of the cumulative 
MRH distributions in average scaling for the two cases, 
A-Frmg- (x) = F RKG (x) - F RMG . (x), where RMG and RMG' de- 
note the stationary and non-stationary ground-states re- 
spectively of the same size L, for L = 128, 256, 512. We 
observe that for all sizes the difference is appreciable, of 
the same order than the difference between the RMG case 
and its Gaussian approximation, and it does not have 
any appreciable decrease with increasing L. We thus con- 
clude that the MRH distribution of RM ground-states is 
sensitive to the boundary conditions. 

In order to understand the origin of the sensitivity to 
boundary conditions we have also compared the MRH 
statistics of Gaussian interfaces described by the sine 
series (defined above) and by the full periodic series of 
Eq. p]). A Kolmogorov-Smirnov statistical test over 10 6 
averaged-scaled numerical Gaussian realizations shows 
that the MRH distribution for the sine series has, for 
( = 1/2, a probability p ~ 0.9 of being the same than 
the MRH distribution of the full-periodic series, while 
for ( = 2/3 the same probability is p ~ 10 -11 . While for 
( = 1/2 these results can be simply related to the fact 
that the signals are Markovian (and thus quickly loose 
"memory" of the border), for £ = 2/3 they show that 
the dependence on boundary conditions is indeed closely 
related to the anomalous geometry, rather than to its 
particular microscopic origin or to the presence of non- 
Gaussian corrections. This extends the results obtained 
for the width distribution of Gaussian signals [16] to the 
case of the MRH. 
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FIG. 10. Differences AF ma ,(x) = -F RM a(a;) — Fbxb>(x) of 
cumulative distributions between the two ways of select- 
ing the optimal path in a disordered medium, for lengths 
L = 128, 256, 512. The cumulative distributions are com- 
puted with average scaling, i.e. with x = h m /h m . 
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V. CONCLUSIONS 

We have studied the maximal relative height statis- 
tics of self-affine one-dimensional elastic interfaces in a 
random potential. We have analyzed, by exact numer- 
ical sampling methods, the ground state configuration 
and the critical configuration at depinning, both for the 
experimentally relevant random-manifold and random- 
periodic universality classes. We have also analyzed the 
MRH distribution of self-affine signals generated by inde- 
pendent Gaussian Fourier modes and obtained an exact 
analytical description for the right tails using Pickands' 
theorem. We found that in general, the independent 
Gaussian modes model provides a good approximation 
for predicting the MRH statistics, as it was previously 
found for the width distribution at depinning. This re- 
sult is a priori not obvious as the MRH is an extreme 
observable (occurring at a single point in space for each 
realization) while the width is a spatially averaged one. 

The comparison of the different MRH distributions 
with their approximation using Gaussian independent 
Fourier modes shows small differences that quickly satu- 
rate with system size in the case of depinning, in agree- 
ment with the results obtained for the width distribution 
of these configurations [HI [T3] : this confirms the predic- 
tions of non-Gaussian corrections at depinning. It would 
be very interesting to characterize more quantitatively 
these non-Gaussian corrections that we have observed for 
the MRH at the depinning threshold using the tools of 
Functional Renormalization Group [14] . A preliminary 
analysis of this question indicates that this is not a sim- 
ple extension of the previous works on the distribution 
of the width [TH Q2] : indeed the extension of the pertur- 
bation theory proposed in Refs. [TU 02] to the computa- 
tion of the MRH distribution involves, at lowest order, a 
constrained propagator (with an absorbing boundary at 
h = 0) for which the Wick's theorem does not hold. 

On the other hand, for ground-state configurations 
we also find small differences but they display a slower 
decrease with system size. For the random-manifold 
ground-states a finite difference is found in the large size 
limit, implying non-Gaussian corrections. This is qualita- 
tively consistent with analytical results obtained for other 
physical observables for the directed polymer in a random 
medium. For instance, the energy of the optimal polymer 
which is described by one of the Tracy- Widom distribu- 
tions (associated to the Gaussian Unitary Ensemble or to 
the Gaussian Orthogonal Ensemble of random matrices 
depending on the geometry of the problem) , which thus 
shows strong deviations from Gaussian fluctuations |38j . 
The random-periodic ground states on the other hand 
have a MRH distribution that is indistinguishable from 
the Airy distribution for the largest system sizes, as fol- 
lows from studying its moments and by a Kolmogorov- 
Smirnov statistical test. This further confirms the uni- 
versality of the Airy distribution for periodic signals with 
£ = 1/2, as shown for a family of non-disordered ther- 
mally equilibrated one-dimensional solid-on-solid mod- 



els [BJ. Our results are also consistent with previous re- 
sults obtained for the related Solid-on-Solid model on a 
disordered substrate in 1 + 1 dimensions (which belongs 
to the same universality class as RPG) [27J . For this 
model with one free end, the ground state configuration 
can be constructed iteratively so that the height field be- 
haves, on large scale, as a random walk with ( = 1/2. 
Notice however that at variance with Ref. [37J we have 
considered here periodic boundary conditions. 

Finally, we have shown that MRH distributions are 
in general sensitive to boundary conditions. This might 
be important for experiments, where particular bound- 
ary conditions are difficult to realize or are not precisely 
known. In this respect it would be interesting to per- 
form an analysis of the MRH statistics using the "window 
boundary conditions" or "segment statistics" proposed 
and applied to the width statistics in Ref. [16] . 

Our results show that the MRH statistics provide a 
valuable tool to study experimental images of self-affine 
interfaces, such as magnetic domain walls or contact lines 
in partial wetting. In particular, it allows to infer infor- 
mation about the mechanisms behind the universal self- 
affine geometry, such as the disorder-induced coupling 
between Fourier modes, or more generally the one in- 
duced by non-linear interaction terms. It also allows to 
infer information about boundary conditions. Our re- 
sults might be used as a guide for analytical approaches 
predicting the statistics of extreme geometrical observ- 
ables. 
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Appendix A: Pickands' theorem and the right tail of 
the distribution of the MRH for a Gaussian path 

Let X(t), t £ [0, T] be a continuous centered Gaussian 
process with covariance function r(t) — (X(s)X(s + t)) 
which satisfies 

1. r(t) < 1, forte [0,T] , 

2. r(t) = 1 - C|*|^ + o(|*H as t ->• , 

where T > 0, n £ (0, 2] and C > are constants. Let us 
define 



X max = max X(t) 

0<t<T 



(Al) 



Pickands' results concerns the right tail of the distribu- 
tion of A max , Pr(A" max < [3T]: 



where $(£) = -4= J. 00 e"" 2 / 2 and _ff M > 0, the so-called 
Pickands' constant, is given by 

H u = lim T _1 / e s Pr f max Y(t) > s) ds . 



F(t)=v^ J B„ /2 (t)-t' i 



(A3) 



where B^^^t) is the fractional Brownian motion with 
Hurst exponent jjl/2, i.e. the Gaussian process charac- 
terized by the two-point correlation function: 

(B, /2 (t)B, /2 (s)) = \ (f + a" - |t - s\n . (A4) 

No explicit formula exist for H^ except for /i = 1 and 
|U = 2, for which H i = 1 and _ff 2 = l/v 7 ^- The result 
above ( A2 1 means that the pdf of X max behaves, for large 



argument as 
d 



je Pr(A max < ~ TC'^^e^e'^ . (A5) 
at, V27T 



Here we will apply this result ( A5 ) to derive the asymp- 
totic behavior of h m = maxo< x <L h(x) where h{x) is 
given in Eq. (fTl) and distributed according a Gaussian 
probability measure as in Eq. (12). We thus first com- 



pute the two-point correlation function as 



(h{x)h{x')) = 



•IL 



Q-l °° 



(2tt 



2irn 



,,, Z^ 71 a C0S i^^-^')), (A6) 



which is obviously stationary. One has in particular 
wf = (h(x)h(x)) — 2 (27r)a C( Q; ); independently of x. If 
one defines 



/i(x) 



fo(x) 

W L 



(A7) 



then the two-point correlator is also stationary r(x) = 
(h(y)h(y + x)) and periodic r(0) = r(L). In addition one 
has from Eq. (A6| that it satisfies r(x) < 1, Vx £ [0,L]. 
To apply Pickands' theorem we need to analyze the small 
x behavior of r(x). To this purpose, it is useful to write 
r(x) as 



r(x) 



2C(«) 



Li 



i(a,< 



fli-KX J L 



Li a,e~ 2 " x/L 



where Li(a, x) is the polylogarithm function [3 
ticular, for a = 2n one has [3"2"| 

r(x) = (-l) n+1 — \4— — B 2n (x) 

w v ; 2(2n)!C(2n) v ; 

B 2n {x) 

~ B 2n (0) ' 



(A8) 
In par- 



(A9) 



where B 2n {x) is the Bernoulli polynomial of degree 2n 
and B 2n (0) = b 2n is a Bernoulli number, which satisfies 
b 2 „ = 2(-l)™- 1 C(2n)(2n)!/(27r) 2n . Hence one has for 
instance, for a = 2 

r(x) = 1 - 6x + 6x 2 , (A10) 

so that, in that case, \x = 1 while for a = 4 one has 

r(x) = 1 - 30x 2 + 60x 3 - 30x 4 , (All) 

so that /i = 2 in that case. The small argument behavior 
of r(x) can be obtained for any a from a careful asymp- 
totic analysis of Eq. ( A8 1 which yields 



r(x) 



2tt 



|r(l-a)sir. 

C(«) 



M, 



1 < a < 3 . 



|7^yX 2 logx , a = 3 



(A12) 



l-2^ 2 ^x 2 ,a>3 



It is then straightforward to use Pickands' theorem (A5) 



to obtain the results given in the text in Eqs ( 19 20 22 



2l|). Note that in the general case where 1 < a < 3, one 



obtains the amplitude D as 
H a -i ( 2tt 



D 



2tt VC(«) 



ri 



a)sin(^) 






(A13) 



where H a _i is the Pickands' constant given in Eq. (A3 1 
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